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ABSTRACT 

The cumulant correlators, C pq , are statistical quantities that generalise the better- 
known S p parameters; the former are obtained from the two-point probability distri- 
bution function of the density fluctuations while the latter describe only the one-point 
distribution. If galaxy clustering develops from Gaussian initial fluctuations and a 
small-angle approximation is adopted, standard perturbative methods suggest a par- 
ticular hierarchical relationship of the C pq for projected clustering data, such as the 
APM survey. We establish the usefulness of the two-point cumulants for describing 
hierarchical clustering by comparing such calculations against available measurements 
from projected catalogs, finding very good agreement. We extend the idea of cumulant 
correlators to multi-point generalised cumulant correlators (related to the higher-order 
correlation functions). We extend previous studies in the highly nonlinear regime to 
express the generalised cumulant correlators in terms of the underlying "tree ampli- 
tudes" of hierarchical scaling models. Such considerations lead to a technique for de- 
termining these hierarchical amplitudes, to arbitrary order, from galaxy catalogs and 
numerical simulations. Knowledge of these amplitudes yields important clues about 
the phenomenology of gravitational clustering. For instance, we show that three-point 
cumulant correlator can be used to separate the tree amplitudes up to sixth order. We 
also combine the particular hierarchical ansatz of Bernardeau & Schaeffer (1992) with 
extended and hyperextended perturbation theory to infer values of the tree amplitudes 
in the highly nonlinear regime. 
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1 INTRODUCTION 

The Newtonian gravitational force is scale-free. This leads one directly to the suspicion that, at least in the strong clustering 
limit where memory of any specific length scale set by the initial conditions has been obliterated, the complex pattern of 
gravity-driven clustering should display relatively simple scaling properties. One way such scaling should manifest itself is in 
the behaviour of the hierarchy of A-point correlation functions (e.g. Peebles 1980). This has led to the formulation of the 
so-called hierarchical ansatz for gravitational clustering which, in its most general form, can be written , 
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(e.g. Balian & Schaeffer 1989), where £jv is the Appoint correlation function for particles located at positions ri ...tn, and 
7 is the negative slope of two-point correlation function £(ri,rj). A more restrictive, but still quite general, algebraic form is 
often used in the literature, in which the correlation function is constructed from linear superposition of all possible topologies 
of "tree" diagrams connecting N points with N — 1 edges. This model permits the association of a number with each distinct 
tree topology (the tree amplitude) with the amplitudes themselves left arbitrary (Balian & Schaeffer 1989). In other words, 
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where the last term is the two-point correlation function for each pair of particles It is important to note that this 

model is not the only possible way to formulate hierarchical clustering. In general, the amplitudes Qjv, a could depend on 
the geometry of a specific configuration of N points, and not just its topology as is assumed here. However, there are some 
indications, at least at the level of three points, that the amplitudes do become independent of configuration in the highly 
non- linear regime (Scoccimarro et al. 1998). There are some indications that they become independent of level of non-linearity 
too at highly non-linear regime (Munshi et al. 1999; see Colombi et al. 1996 for a different view-point). 

Many different models of hierarchical scaling, largely consisting of particular choices for the Qn,o, can be found in the 
literature. In principle, it should be possible to derive these amplitudes from the underlying gravitational physics encoded in 
the hierarchy of BBGKY equations that describe the dynamical origin of gravity-driven clustering (Davis & Peebles 1977; 
Fry 1984; Hamilton 1988). However, no general solution to these equations exists, so these various approximate models are 
employed in order to get some insight into the physics of gravitational clustering in highly nonlinear regime and to aid 
interpretation of empirical galaxy clustering data. 

Despite the increasing availability of large samples of galaxy redshifts, the most precise statistical analyses of galaxy 
clustering are still performed using projected (angular) samples. The extremely large size of projected galaxy surveys, together 
with the relatively easy-to-handle effect of projection, makes such samples ideal for the extraction of descriptors such as 
higher-order cumulants. The APM survey (Maddox et al. 1990; Maddox, Sutherland & Efstathiou 1990; Maddox, Efstathiou 
& Sutherland 1996), which contains more than 1.3 million objects, is at present the largest angular catalogue available. Such 
is the amount of useful information in this sample that Gaztanaga (1994) was able to extract estimates of the so-called 
S p parameters (see below), for p < 9 from it. Such an analysis is at present beyond the scope of galaxy redshift surveys 
(Bouchet et al. 1993), at least until the Sloan Digital Sky Survey and Anglo- Australian 2DF survey are completed. Projected 
catalogs are also free from redshift-space distortions, a fact which can be exploited to yield accurate estimates of real-space 
clustering phenomena in a relatively unambiguous way. The obvious difficulty associated with projected catalogs is related to 
the fact that any projected angular scale corresponds to a superposition of physical length scales which have which represent 
different levels of nonlinearity. Limber (1954) pioneered the study of relating angular correlation functions to their real-space 
counterparts. These inversion techniques and their small angle approximations are now known to be very efficient for the 
extraction of information from projected catalogs (e.g. Peebles 1980; Baugh & Efstathiou 1993). 

Many statistical studies of large-scale structure focus (in one way or another) on the determination of one-point cumulants 
for cells of a particular size, which are essentially volume averages of correlation functions over the cell. The second-order 
cumulant, for example, is the variance of cell-counts, and which is the average of the two-point correlation function over the 
cell. The standard S p parameters (see below) are based on the cumulants of the one-point probability distribution of density 
fluctuations, fi(5), and can be expressed in terms of volume averages of the higher-order correlation functions in an analogous 
fashion. Here and throughout the paper <5(x) denotes the dimensionless density contrast, defined by 

= £(x)- £ o j 

Po 

where po is the mean density of matter. 

Szapudi & Szalay (1997) recently proposed a new class of statistical measures, known as the cumulant correlators, which 
generalises the set of S p parameters mentioned above. The cumulant correlators C pq are analogous quantities obtained from 
the two-point (bivariate) distribution f 2(61, 62). They form a symmetric matrix, C pq , and estimation of these quantities from 
data is a very similar task to the extraction of the standard two-point correlation function. Furthermore, using the factorial 
moments of the discrete count probability density function makes it possible to subtract the effect of Poisson shot noise very 
efficiently. Szapudi & Szalay (1997) extracted estimates of these quantities for p + q < 5 from the APM galaxy survey. 

Theoretical work on the cumulant correlators has been performed by Bernardeau (1995), who obtained predictions for 
the C pq in three-dimensional real space for clustering developing from Gaussian initial conditions and in the limit of large 
spatial separations. Interestingly, a relatively simple scaling hierarchy is expected to develop under these conditions, identical 
to the scaling of their one-point counterparts the S p parameters, potentially furnishing a simple and powerful diagnostic of 
the idea of gravitational instability and, perhaps, of the extent to which galaxy bias affects clustering statistics. 

In this paper we study the generalisation of the cumulant correlators to multi-points and for arbitrary order. These 
statistics are constructed from the quantities (<5 p (xi)<5 9 (x2) . . . <5 s (x;)), but while the cumulant correlators have I = 2 with p 
and q arbitrary, the multipoint correlators have no restriction on I or any of the integers p, q, . . . s. We study the dependence 
of these quantities on the underlying hierarchical tree amplitudes, with a view to testing the different possible hierarchical 
models using these statistics. We also present the results of a calculation of the quasilinear behaviour of the cumulant 
correlators in projection, including the effects of observational selection. In Section 2 we outline our main theoretical results 
in highly nonlinear regime; Section 3 is devoted to the computation of two-point cumulant correlators in projected catalogues 
in quasi-linear regime. Theoretical results are tested against measurements from the APM survey in Section 4. 



2 GENERALIZED CUMULANT CORRELATORS IN THE NONLINEAR REGIME 

Such are the importance of hierarchical models in the understanding of gravitational clustering that it is important to learn 
which of choice of tree amplitudes in the model (2) corresponds with reality. This is best tackled by attempting to extract these 
amplitudes from numerical simulation data. However, the direct determination of higher-order correlations is complicated by 



their possible dependence on shape as well as topology. Most numerical studies in this direction are therefore have concentrated 
on an indirect determination these quantities, using one-point cumulants or two-point cumulant correlators. It is known that 
cumulant correlators can be used to compute the amplitudes associated with two different types of topologies - the so-called 
"snake" and "star" topologies - at fourth order (Szapudi & Szalay 1997). However, at higher orders than this the number of 
different topologies rapidly increases and two-point cumulant correlators are no longer sufficient to effect a unique separation 
of all tree-amplitudes. The only possible way we can remedy the situation is by moving to the generalised cumulant correlators 
mentioned above. These are the natural generalisation of their two-point counterparts (Munshi et al. 1998a) and, as we discuss 
elsewhere, they are also related in a very interesting way to the statistics of collapsed objects, in much the same way as the 
two-point cumulant correlators are related to the natural bias associated with treating collapsed objects rather than the 
general distribution of mass (Bernardeau 1992; Munshi et al. 1998a). 

The terminology in this field is a little confusing, so it is worth beginning with a brief explanation of notation in order 
to relate this work with the existing literature. For three points there is only one possible tree, and the unique three-point 
amplitude is generally denoted Q. As we mentioned in passing, there are two possible topologies for connecting four points 
(respectively called the "snake" and the "star"). The amplitudes of these tree configurations are usually denoted R a and Rb- 
Following in this vein we introduce, for the distinct topologies at fifth order, the notations S a , S c and Sb for the "snake", 
"star" and "hybrid" cases. At sixth order the number of different topologies is six (Fry 1984) and one can similarly use T a , 
Tb, T c , Td, T e and T/ to denote the amplitudes associated with them. 

The order M of the multi-point cumulant correlator (c5 p (xi)<5 9 (x 2 ) . . . 5 s (x;)) is defined to be equal to p + q + . . . + s 
and such a quantity of order M takes contributions from M-point correlation function. Cumulant correlators of a particular 
order will depend only on hierarchical amplitudes associated with topologies contributing to the same order, even though the 
number of spatial points (I) on which they depend may vary. 

For example, the sixth-order two-point and three-point cumulant correlators will both depend on tree amplitudes as- 
sociated with six-point correlation functions. For example, there are three different but "degenerate" three-point cumulant 
correlators of sixth order, a fact which stems from the number of different ways the integer six can be decomposed in three 
integers, i.e. 4+1 + 1 = 3 + 2 + 1 = 2 + 2 + 2 = 6. All of these three-point cumulant correlators in general can depend on all the 
six different tree-amplitudes associated with sixth order correlation function. Hence measurement of these three three-point 
cumulant correlators of sixth order will provide us three equations connecting the six tree amplitudes. The system of three 
equations with six variables is consequently indeterminate, and can only be solved by including three more equations connect- 
ing these amplitudes. But if we compliment these equations with information provided by the two-point cumulant correlators 
of sixth order we get three more equations, corresponding to the fact that 1 + 5 = 2 + 4 = 3 + 3 = 6 and quantities also depend 
on exactly the same tree amplitudes. The resulting six equations in six variables will provide us with a fully-determined set 
of equations which can be solved to yield the amplitudes associated with all topologies of sixth order. 

This example can be generalised. Counting the different number of degenerate multi-point cumulant correlators, for a 
given order and for a specific set of points, one can check that three-point cumulant correlators can in principle be used to 
determine all tree amplitudes associated with correlation functions up to and including sixth order. We have tabulated the 
different multi-point cumulants of at a given order to show which can be used to separate out the hierarchical amplitudes 
associated with correlation functions through equation (2) . 



2.1 Leading-order tree terms 

If we express two degenerate two-point cumulant correlators of fourth order in the limit of large separations (or, in other 
words, the variance evaluated at each point, related to £ 2 , is much bigger than the covariance between the points, expressed 
as £ a b) we obtain following pair of equations: 

(«5 3 (xi)<5(x 2 )> = (3R b + 6R a )U& 2 ; (4) 
(«5 2 (xi)(5 2 (x 2 )) =4R b U& 2 (5) 

As explained earlier these two-point cumulant correlators depend only on the two amplitudes associated with the four-point 
correlation function, i.e. R a and Rb- They can solved simultaneously to compute values of these quantities. 

At fifth order there are also two degenerate two-point cumulant correlators which depend on three tree-amplitudes: 

(5 4 ( Xl )5(x 2 )) = (4S C + 36S 6 + 24S a )U& 3 (6) 

(5 3 (x 1 )5 2 (x 2 )) = (6S b + 12S.)U6 3 . (7) 

To determine the amplitudes S a , Sb and S c the two-point cumulant correlators would be insufficient, so we have to move to 
the three-point quantities instead. These provide us another pair of equations with same tree-amplitudes: 

(5 2 (xi)<5(x 2 )<5 2 (x 3 )) = {4SaU&c + 4S 6 £ a6 £ ac + 4S' f) £ ac £ i)c )f 2 2 (8) 
(5 3 (xi)5(x 2 )5(x 3 )> = {{6S a + 3S b )(Utibc + £«<&<:) + (3S C + lSS b + 6S a )£ a6 £ ac )6 2 (9) 

Any three of the equations (6) to (9) can be used to determine the tree-amplitudes S a ,Sb and S c - The other equation will 
be consistent but redundant, and can be used as a check. Of course, the one-point cumulant of same order always provide an 
additional constraint at every order. 



There are three degenerate two-point cumulant correlators of sixth order and there are 6 distinct topologies: 

(<5 4 (xi)5 2 (x 2 )) = (48T a + 48T e + 24T 6 + 8T d )£ ab £ 2 4 (10) 

<<5 3 (xi)<5 3 (x 2 )) = (36T 6 + 36T a + 9T C )U& (H) 

(5 5 (xi)<5(x 2 )> = (120T a + 180r 6 + 180T e + 80T d + 60T C + 5T f )£ ab £ 2 4 (12) 

Combining these three equations with expressions for the three-point cumulant correlators at sixth order we can have another 
set of three equations which can be used to used to determine all the relevant tree-amplitudes. 

(<5 2 (x 1 )«5 2 (x 2 )5 2 (x 3 )) = (8T a + 8T e )(U6c + + UZbc)b 3 (13) 

(<5(xi)5 2 (x 2 )5 3 (x 3 )} = ((18T 6 + 12T a + 6T c )£ a6 &c + (12T„ + 24T b + 12T e + 6T d )£ a6 £ ac + (6T 6 + 12T a )£ oc £ a6 )£2 3 (14) 

(<5(xi)5(x 2 )<5 4 (x 3 )) = (24T a + 24T e + 12T 6 + 4T d )(£ a6 £ 6c + UZac) + (4T, + 24T e + 36T C + 48T d + 24T a + 120T 6 )£ OC 6 C )£ 2 3 (15) 

This method can be extended to higher-order correlation functions and, combined with a judicious use of factorial correlators 
or factorial moments described in Munshi et al. (1998b), it provides the simplest way to determine the hierarchical amplitudes. 
It is however to be noted that we have neglected the small correction factors pointed out by Boschan et. al. (1994). 

2.2 Generalised Cumulant Correlators in the Bernardeau— Schaeffer Ansatz 

As explained earlier, although all hierarchical models (2) agree on the basic tree structure underpinning the correlation 
hierarchy describing the matter distribution, specific models differ from each other in the assumptions they make regarding 
the amplitudes associated with tree-topologies. One such model, which has been studied in great detail, was proposed by 
Bernardeau & Schaeffer (1992). They assumed that each vertex (or node) in the tree representation of correlation functions 
(2) carries an weight which depends only on the order of the vertex and on nothing else. Vertices of the same order appearing 
in different trees will carry exactly same weight. The amplitude of the whole tree is consequently always equal to the product 
of the weights assigned to each constituent node. Even with this rather simple assumption this model can predict many 
interesting features associated with collapsed objects, some of which have already been tested successfully against numerical 
simulation (Munshi et al. 1998b). More sophisticated models have also been constructed by assuming a specific form for the 
generating functions of these vertex amplitudes. Such models can make further predictions about many body statistics of 
collapsed objects. 

Assuming that the tree amplitudes can be factorized in this way, it is possible to decompose the cumulant correlators of 
various orders in the following way (Bernardeau 1992; Munshi et. al. 1998a): for the two-point function of arbitrary order 

<<5 p (x a )5 9 (x ft )} = [C vl UC q i\e +q - 2 ; (16) 
for three points, 

(<5 P (x a )(5' J (x i ,)(5 r (x c )) = [Cql^abCpll^acC'rl + Cpl^abCqll^bcCrl + Cpl^acC'rll^bcCplj^ P + q+r 3 '; (17) 

for four points 

(S p (x a )8 q (x b )5 r (x c )5 s (x d )} = [C q x^ ab CpXYLiadCsiia C C r i + (cyclic permutations) 

+C q i£abC p ii£acC r ii£cdCsi + (cyclic permutation)] £P+i+ r + s - 4 - (ig) 

and for five points 

((5 p (x )J 9 (xi,)(5*"(x c )(5 s (x d )5*(x e )} = [Cpim^abCqi^acCri^aeCti + (cyclic permutation) 

+ C p iii£ a dC q ii£ bc C r i£adCti£aeC 3 i + (cyclic permutation) 

+ CpitabCgiitbcCriitcdCsiitdeCn + (cyclic permutation)] C P+9+r+s+t_5 (19) 
The vertex functions C p i...i are simplest to deal with when expressed in terms of their generating functions, e.g. 

oc 

P =i 

CO 

P =i 

and so on. With these definitions it is possible to relate the generating functions ty p (ti, . . . ,t p ) of the cumulant correlators 
with generating functions of the vertex weights (see Munshi et al. 1998a for more details): 

* (2) (ti,t 2 ) = m(-ti)UM-t2) (22) 
* (3) (ti,t 2 ,t 3 ) = Mi(-£i)£a&£t2(-t2)£ac£*i(-t 3 ) + ... (cyclic permutations) (23) 



* (4) (ti,t 2 ,i3,i4) = Mi(-*i)£a&M3(-£2)£ac£*i(-t3)£aci£*i(-i4) + • • ■ (cyclic permutations) 

+Vi{-ti)iab^{-te)£acli2{-tz)( > adlJ-i{-U) + ■ ■ ■ (cyclic permutations) (24) 

Specific forms for fi n (t) have been obtained both by diagrammatic techniques and direct analytical computation in Munshi 
et al. (1998a). Direct evaluation of the multi-point cumulant correlators provides us with a way to check such predictions to 
arbitrary order. For example, according to the model of Bernardeau & Schaeffcr (1992), the topological amplitudes at higher 
order are always related with the lower order amplitudes in a specific way: at every order this particular model introduces a 
new "star" diagram, and all other diagrams at every order can be expressed in terms of lower order amplitudes in a simple 
fashion. For example, in this model, R a (4th order) is generated by Q — Q3 (third): 

Ra = Q 2 , (25) 
while R b is a "new" amplitude at fourth order, which cannot be related to Q. At fifth order, 

S a = Q 3 ,S b = QR b (26) 
but S c is not constrained by Q or R a or R b . Likewise, 

T a = Q 4 , T b = Q 2 R b , T c = Rl T d = QS C , T e = Q 2 R b , (27) 

with Tf free at this level. Tests of these predictions will become possible as larger simulation boxes become available. These 
will undoubtedly lead us to better understanding of how this model describes the nonlinear regime of gravitational clustering. 



2.3 Generalised Cumulant Correlators in the Szapudi— Szalay Ansatz 

An analysis similar to the preceding one has been carried out by Szapudi & Szalay (1993a), using a different ansatz for the 
tree-level amplitudes to that of Bernardeau & Schaeffer (1992). Szapudi & Szalay (1993b) subsequently also looked at the 
predicted statistics of over-dense cells in this light. In their analysis they assume that amplitudes associated with different 
topologies but with same number of vertices are always equal. For example this means that, at fourth order, all tree topologies 
will have an amplitude Q4 which is a weighted average of different topologies of fourth order. Since there are 16 distinct 
configurations of 4 points, 12 of which are "snake" and four of which are "star", this means I6Q4 = 12R a + 4R b . Similarly, 
at fifth order, 125<5s = 60S a + 605*6 + 5S C and in general Qn,c = Qn- The corresponding formulae are: 

{5 p (* a )8 q {* b ))=p p - 1 q' l - 1 Q J>+q Ue +q ~ 2 (28) 

(«F(x a )<5*(x 6 )«f (X C )) = Y" V^Q^+r [ P Uiac + qUZbc + V^bc] f +9+I - 3 (29) 

(5 p (x a )<5 9 (x(,)<r(x c )<5 s (x d )) = p p ~ 1 q q ~ 1 r r ~ 1 s s ~ 1 Q p+q + r + s [p 2 ^ab^ ac ^ad + (cyclic permutations) 

+ priabia^cd + (cyclic permutations)] I* +i+ r + s ~ 4 (30) 

{8 p (x a )S q (x b )8 r (x c )8 s (x d )S t (x e )} = p p ~ 1 q q ~ 1 r r ~ 1 s s ~ 1 t t ~ 1 Q p+q+r + s +t [p 3 £ a b£ac(, a d£,ae + (cyclic permutations) 

+ P 2 q(,ab£,bc(,ad£,ae + (cyclic permutations) 

+ qrs^ ab ^ bc ^ cd ^ de + (cyclic permutations)] ^P+i+ r + s + t ~ 5 (31) 

These results again furnish a simple test of this model. One simply needs to evaluate Q p +...+t for different values of the 
particular p, q . . . ,t and check if the parameter Q p +...+t remains invariant when the individual p, q . . .t are changed in such 
a way that p + . . . + t remains constant. 

Using two-point cumulant corrrelators, Munshi & Melott (1998) have recently tested predictions of these particular 
versions of the generic hierarchical ansatz against numerical simulations up-to fourth order. More detail studies using three- 
point factorial correlators up to sixth order will be presented elsewhere. 



2.4 Hierarchical Amplitudes from Extended and Hyper-extended Perturbation Theory 

Extended perturbation theory and hyper-extended perturbation theory have been suggested by Colombi et al. (1996) and 
Scoccimarro & Frieman (1998) respectively in attempts to obtain more accurate analytic predictions of the 5jv parameters 
in the highly non-linear regime. Such extensions however are, however, silent concerning the specific contribution to Sn from 
distinct topologies. We can improve upon this situation, by combining predictions of the hierarchical model of Bernardeau & 
Schaeffer (1992) for the highly non-linear regime with these predictions using (hyper)extended perturbation techniques. 

It turns out that this procedure always predicts that "star" topologies carry more weight than "snake" and "hybrid" 
topologies. We also find that the amplitudes of all topologies increases when the fluctuations have more large-scale power. For 
power spectra with less power on larger scales we find that all topological amplitudes become approximately equal, although 
perturbation theory is known to break down for such spectra. Existing numerical results at fourth order (R a and R b ) roughly 
match with these predictions (Munshi & Melott 1998). However, more tests are needed to check predictions at orders greater 
than four. 
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Figure 1. Predictions for hierarchical amplitudes of fourth order R a , and Rf, and fifth order S a , Sb and Sc are plotted as a function 
lowest order hierarchical amplitude Q. Solid lines are predictions from hyper-extended perturbation theory (Scoccimarro & Frieman 
1998), dotted lines represent extended perturbation theory (Colombi et al. 1996) and dashed lines represent predictions from specific 
hierarchical model of Bernardeau & Schaeffer (1992). The bottom most curve in each panel represents snake diagrams R a = Q 2 and 
•S a = Q S which are same for all these approximations are same. The triplets of upper curves in left panel represent R b and in right panel 
they represent S c and Si, respectively from top to bottom. 



It may not be possible to extend tree-level perturbation techniques straightforwardly by changing spectral index n to 
effective spectral index n ft in order to compute multi-point cumulant correlators. In principle, although these quantities can 
be computed for arbitrary points in the quasi-linear regime they will in general imply that hierarchical amplitudes are depend 
on shape parameters while, in the highly nonlinear regime hierarchical amplitudes are thought to be independent of shape. At 
least at the level of three-point this seems to have been already confirmed by some numerical experiments (Scoccimarro et al. 
1998). Furthermore, it has been pointed out (Bernardeau 1996) that, in the quasilinear regime, tree-level perturbation theory 
for one-point cumulant of the smoothed density field 53, implies 53 = 3Q but for C21 is not equal to 2Q as one would have 
expected on the basis of the hierarchical ansatz (which applies to the unsmoothed tree- level perturbation theory). Clearly, 
one has to be careful about the effect of smoothing on the cumulant correlators. When combining the hierarchical ansatz 
with extended or hyperextended perturbation theory we have to keep these points in mind. It is nevertheless interesting to 
combine these two kinds of calculation in this way, as there are otherwise no analytical predictions for these quantities that 
can be tested against numerical simulations. 

The results of this and the previous calculations are shown in Figure 1. 



3 QUASI-LINEAR CUMULANT CORRELATORS FROM PROJECTED CATALOGS 

In the quasi-linear regime, two-point cumulant correlators have already been studied in great detail using perturbation theory. 
The extension of such studies to multiple points, as is required for the generalisations we discuss in this paper, are difficult 
owing to the complicated shape dependence of the hierarchical amplitudes. For this Section, therefore, we shall focus on the 
(original) two-point version of these statistics. 

Estimates of the cumulant correlators were extracted from the APM catalogue by Szapudi & Szalay (1997) and a related 
theoretical computation of projected cumulant correlators was done by Bernardeau et al. (1997). Let us now concentrate on 
the APM results in order to illustrate the excellent match between theory and observation. 

The usual notation with which the two-point cumulant correlators are expressed is in terms of C pq , axe defined as 
(Bernardeau 1995; Szapudi & Szalay 1997) 

(^(xi)^(x 2 )') c = C P9 (5 2 (x)>(5(x 1 )«5(x 2 )) p+9 - 2 . (32) 



An alternative parameterisation, which we shall in fact use, is in terms of the quantities 



We assume the initial density contrast <5(x) (3) to be Gaussian which means that it can be decomposed into independent 
Fourier modes <5(k) 

<5(x) = y*d 3 k(5(k)exp(ik.x), (34) 
which completely characterized by power spectra P(k) 

(<5(k)<5(k')> = 6 D (k + k')P{k). (35) 

In such a field, each Fourier mode is statistically independent of the others. 

By analogy, for projected catalogs, one define angular cumulant correlators to be 

M 7l )a( 72 )) c = Cp9 <a 2 ( 7 )>M 7l M 72 )) p + 9 - 2 , (36) 

where 7 i and 7 2 are unit vectors defining positions in the celestial sphere. A quantity q can be defined with respect to c in 
the same way as analogy with Q is defined with respect to C in equation (33). The quantity c pq is the projected cumulant 
correlator and a(j) is projected number density of objects: 

(7(7) = no / r 2 drF(r)^^-, (37) 
Jo P° 

where F(r) is the selection function of the catalog and po the mean density of the universe, noF(r) is mean density of 
observable objects at a radial distance r. The projected density field on the sky will usually be smoothed with a smoothing 
window W 70 , the result of which can be written as 

M 7 ) = n / ^7o(7i)rf 2 7i J r 2 drF{r)^-. (38) 

The function W / ( 7 o) is unity if the angle between 7 and a given direction 7o is less than 8o and zero otherwise. The projected 
one-point cumulants 

Sv = K(7)>/(^(7)) p - 1 (39) 

are normalised moments of the smoothed one-point projected probability density functions, and coincide with the cumulant 
correlators in the limit xi — > x 2 . 

The most important assumption that we will be using in our derivation of c pq in projection is the fact that both smoothing 
angle and the separation angle are small compared to unity. In general, c pq will contain terms which are of higher order in 
(cr( 7 i)a"( 7 2)}, but we will be neglecting all such terms in our calculations. Detailed perturbative calculations are needed to take 
higher order loop corrections into account. In this paper we however focus on directly comparing these theoretical predictions 
against measurements from the APM (Szapudi & Szalay 1997). One of the advantage of using small angle approximation 
is the separation of effects due to selection function and that of dynamical contributions. Using these approximations, it is 
possible to show that the projected cumulant correlators can be written: 

c pq = R p+q <d pq ; q pq = c pq /q q ' 1 p p ' 1 ; (40) 

where Q pq are the cumulant correlators derived from the simplified dynamics of two-dimensional isotropic collapse in a manner 
we shall shortly describe. We can define the generating function /?(j/i,j/2) for pq as 

00 p q 

/%!,»)= E e -||- ( 41 ) 

P=1,<3=1 

Since we are interested only in terms linear in (<5( 7 i)<5( 7 2)) we can write /3 (3/1,3/2) = (3(yi)f3(y2) (Bernardeau 95), where 

/%i)=f> P if, (42) 
P =i 

which means that, in linear order, Q pq — Q p iQ q i. Bernardeau (1992) showed that the generating function, 



P =i 

of the tree amplitudes v p that appear in a perturbative expansion of the density distribution satisfy dynamical equations 
governing collapse of spherical over density, with r playing the role of time (Bernardeau 1992; Munshi et al. 1994). It can be 
shown that generating function for cumulant correlators (3{y) can be related to G(r) by the following relation (Bernardeau 
1995): 

/%) = -2/«; (44) 




r = -y-^. (45) 

These relations hold only for an unsmoothed initially Gaussian field. In the case of angular smoothing and power law initial 
condition P{k) = Ak n one can write 

/3(j/)=r(y)[l + G(r(y))n (n+2)/2 . (46) 

It is easy to see that although the dynamics is determined by two dimensional spherical collapse, the spectral index n is still 
that of three dimensional power spectrum P(k). 

Expanding P(y) in Taylor series we obtain the dynamical contribution to the cumulant correlators (Bernardeau et al. 
1997; Munshi & Melott 1998) 

21 = y-i(n + 2) (47) 

8« = ^-^(» + 2) + §(n + 2)» (48) 

The factors R p + q in (40) depend on the spectral index and specific form of selection function. They were determined by 
Gaztanaga (1994) 



mr 2 (3) M p (3p-(p-l )(n + 3)) 
Mf-^S-n) 

where 



^ a fP— 1 /o \ V / 



M p (a) = / dr r a ~ 1 F p (r). (50) 



Using the specific form for the selection function F(r) = Kr~ ' 5 exp[— (r/D) 2 ] gives Rz = 1.19, Ra = 1.52, _Rs = 2.00, as 
computed by Gaztanaga (1994) and Bernardeau (1995). The slope of the projected two-point correlation function used for 
comparison is 7 = 1.7. 



4 COMPARING ANALYTICAL CALCULATIONS WITH APM MEASUREMENTS 

Fully nonlinear calculations of q pq using factorial moments were performed by Szapudi & Szalay (1997), who used 0.23° cells 
to construct density maps with a magnitude cuts of bj — 17 to 20. Measurements of q pq were presented up-to fifth order. 
For these magnitude cuts, the angular separation of 1° corresponds approximately to 7/i -1 Mpc, a regime in which where 



perturbative calculations should be valid. We have replotted q pq measured from the APM survey by Szapudi & Szalay (1997) 
as a function of separation angle 0° . The smoothing angle remains fixed for all separation angles. 

The theoretical predictions discussed above are plotted as straight lines and show reasonable agreement with measured 
values. Our results are valid in the limit of large separations and in the limit of small angles. These approximations need to 
be tested with more detailed perturbative calculations and Monte-Carlo simulations. It is clear that dominant contribution 
will come from the tree level approximation used by us. Simultaneous use of these two approximations may appear self 
contradictory, but it has been shown by direct Monte-Carlo integration of ^(ri, r2, ^3) in real space that the value of C21 
converges to its value in the large separation limit, even when the separation of the two cells is only equal to their diameter 
(Bernardeau 1995). It has also been shown (Gaztanaga & Bernardeau 1997), that the small-angle approximation is very good 
for smoothing scales less than 1° and given the errors associated with measurements from galaxy catalogs are large, it provides 
a very good approximation even for separation as large as 5° . On the other hand, it is expected that contributions coming from 
terms of higher order in {a (71) (7(72)) will increase the tree-level values computed here. Finite sample volume corrections tend 
to reduce the computed value of cumulant correlators and enhance the agreement with tree-level perturbative calculations. 

It will be very interesting to include effects of finite volume corrections in measurement of cumulant correlators from 
galaxy surveys and compare with perturbative calculations beyond tree level. However, the important point we wish to stress 
here is that projection effects make it difficult to deduce the values of hierarchical amplitudes, such as R a and Rb, directly. This 
means that a comparison with predictions made in real space and in three dimensions with projected cumulant correlators is 
highly problematic (e.g. Szapudi & Szalay 1997). 

It is also interesting to note that, although the smoothing length involved in this measurement is in the quasi-linear 
regime, the non-linear correction process adopted by Szapudi & Szalay is based on a hierarchical ansatz that applies in the 
highly nonlinear regime. This nevertheless seems to work very well as the values of r a and r b (the projected amplitude of two 
different topologies) derived show a constant amplitude over all angular separations after correction. However, there may be 
an error in their paper. They have derived the relations r a = 522 and r b = 3^31 — 2gi2 for large angular separation 9, as all 
other terms higher order in £;/£ s are negligible in this limit. Taking a closer look in their figure we can see that while the 
measured value of 922 for large angular separations in the lower panel is very close to unity the measured value for 931 in the 
same limit is close to 2.23 (approximately 0.35 in logarithmic units). Using these values of ^31 and q2i we get r a — §22 — 1 
and rb = 2>qsi — 2qn — 3 x 2.23 — 2 x 1 ~ 4.69. However the measured r a remains close to 5.5, and the curve which they label 
rt remains close to unity. This indicates that they have probably inadvertantly interchanged r a and n,. 



5 DISCUSSION 

We have shown how to generalise the concept of a cumulant correlator to an arbitrary number of points, and have deconstructed 
the resulting functions explicitly up to sixth order in large separation limit. These results were obtained without making any 
specific ansatz for amplitudes associated with different tree-topologies. These results, when combined with results of multi- 
point factorial correlators derived in Munshi et al. (1998b), allow us to build a direct determination of the hierarchical 
amplitudes. They can also be used to test the dependence of these parameters on shape parameters in the highly nonlinear 
regime. 

Earlier studies by Szapudi & Szalay (1997) showed that two-point cumulant correlators can separate contributions from 
tree topologies at fourth order. In this paper we have explicitly shown that if we move from two-point to three-point cumulant 
correlators it will be possible to separate all tree topologies contributing to fifth and sixth order correlation functions. For 
higher orders it is necessary to move to four-point cumulant correlators, and so on. 

We have also studied the predictions of more particular versions of the hierarchical ansatz which generally assume some 
specific form for the hierarchical amplitudes. In particular, we have studied an extension of the ansatz given by Bernardeau & 
Schaeffer (1992) which assumes that, in the highly nonlinear regime, a given tree amplitude can be constructed by multiplying 
vertex amplitudes constituting the tree. We have also extended the ansatz by Szapudi & Szalay (1992) which is based on 
replacing each amplitude of the different tree-topologies by an average over different topologies of same order. With our 
analysis it will be possible to test these results against numerical simulations when simulations with larger dynamic range 
become available. This will provide an unique way to study gravitational clustering in the highly nonlinear regime. 

While using our results of multi-point cumulant correlators derived in highly nonlinear regime it should be realized that 
although the decomposition is still valid, the amplitudes depend strongly on shape factors. The shape dependence of lowest 
order tree amplitude Q3 has been studied extensively by Scoccimaro et al. (1998). Similar studies of higher order amplitudes 
will clarify weather these quantities do become independent of shape factors and also more studies are needed to check how 
they depend on initial power spectra. 

We have shown that predictions made for two-point cumulant correlators match very well when compared with mea- 
surement from projected catalogs in the quasilinear regime. Our results in the highly nonlinear regime will also be interesting 
when larger three-dimensional galaxy catalogs are available. It will possible to separate lower order tree topologies using our 
method and to test the validity regime of different nonlinear approximations. 

Hierarchical amplitudes are also related with the statistics of collapsed objects, i.e. how the over-dense cells are distributed 
against background matter distribution (Bernardeau & Schaeffer 1992; Munshi et al. 1998a,b) so study of statistics of collapsed 
objects provide another interesting way to determine the hierarchical amplitudes. The reader is referred to Munshi et al. 
(1998a,b) for further information about this work. 
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Table 1. Multipoint Cumulant Correlators 
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